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Spin-imbalanced ultracold Fermi gases have been widely studied recently as a platform for ex¬ 
ploring the long-sought Fulde-Ferrell-Larkin-Ovchinnikov superfluid phases, but so far conclusive 
evidence has not been found. Here we propose to realize an Fulde-Ferrell (FF) superfluid without 
spin imbalance in a three-dimensional fermionic cold atom optical lattice, where s- and p-orbital 
bands of the lattice are coupled by another weak moving optical lattice. Such coupling leads to a 
spin-independent asymmetric Fermi surface, which, together with the s-wave scattering interaction 
between two spins, yields an FF type of superfluid pairing. Unlike traditional schemes, our proposal 
does not rely on the spin imbalance (or an equivalent Zeeman field) to induce the Fermi surface 
mismatch and provides a completely new route for realizing FF superfluids. 


PACS numbers: 03.75.Ss, 67.85.-d, 74.20.Fg 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, 
characterized by Cooper pairs with finite center-of-mass 
momenta [1, 2], is a central concept for understanding 
many exotic phenomena in different physics branches [3] . 
A crucial ingredient for realizing FFLO states is a large 
Zeeman held that induces a Fermi surface mismatch of 
two paired spins [1, 2]. In recent years, FFLO states 
have been extensively studied in ultracold Fermi gases, 
where the population imbalance between two atomic in¬ 
ternal states (pseudo-spins) serves as an effective Zee- 
man held [4-14] . Despite the intrinsic advantages of cold 
atoms compared to their solid state counterparts, conclu¬ 
sive evidence of FFLO states has not been found yet be¬ 
cause of various obstacles. For instance, a large Zeeman 
held suppresses the superfluid order parameter, leading 
to a very narrow parameter region for FFLO states in 
2D or 3D which can be easily destroyed by thermody¬ 
namic fluctuations [4-6]. In ID, the parameter region for 
FFLO states could be large, but the quantum fluctuation 
is strong [7, 12, 14]. The recently proposed schemes us¬ 
ing spin-orbit coupling and in-plane Zeeman held in a 3D 
Fermi gas may potentially overcome these obstacles [15- 
21] in principle, but they face practical experimental is¬ 
sues such as the large spontaneous photon emission from 
the near-resonant Raman lasers [22-31] and the strong 
three-body loss at Feshbach resonance in the presence of 
spin-orbit coupling [23-26]. 

In this Letter, we propose a new route for realizing 
FF superhuids in ultracold Fermi gases without involv¬ 
ing population imbalance of two spin states that interact 
for generating Cooper pairing. Instead, we induce an 
asymmetric Fermi surface for the generation of FF states 
by other means and the populations of the two spins are 
fully equal. Our main results are the following: 

1) We show that the s- and p x -orbital bands of a 
3D static optical lattice can be coupled using a weak 
ID moving optical lattice along the x direction, which 
can be generated by two counter-propagating lasers 


with the frequency difference matching the s-p x band 
gap. The s- and p^-bands can be denoted as the band 
pseudospin, and the moving lattice induces a band- 
pseudospin-momentum (i.e., spin-orbit) coupling and an 
in-plane Zeeman field, which yield an asymmetric Fermi 
surface along the x direction. The realization of such 
band-pseudospin-momentum coupling may provide a new 
platform for exploring exotic spin-orbit coupling physics. 

2) We show that the asymmetric Fermi surface, to¬ 
gether with the s-wave pairing interaction between two 
equally populated hyperfine spin states, can induce an 
FF type of Cooper pairing within a large parameter re¬ 
gion in the 3D optical lattice, in sharp contrast to the 
narrow parameter region for the spin-imbalanced Fermi 
gas [6, 7]. Because of the 3D nature of the FF superfluids, 
the quantum fluctuations are also suppressed. The gen¬ 
erated FF state is thermodynamically much more stable 
than the spin-imbalanced Fermi gas. Compared to the 
spin-orbit coupled schemes [15-21] that require near res¬ 
onant Raman lasers [23-31], all lasers used here are far- 
detuned, therefore the proposed scheme should work for 
all types of fermionic atoms, including 6 Li [14]. Further¬ 
more, because the hyperfine spins are not coupled with 
the momentum, the s-wave scattering interaction should 
be the same as regular Fermi gases without significant 
three-body loss at Feshbach resonance. These intrinsic 
advantages of our spin-balanced scheme make it experi¬ 
mentally more feasible than the spin-imbalanced schemes 
(with [4-14] or without spin-orbit coupling [15-21]), and 
thus may open a new route for observing FF superfluids. 

Asymmetric Fermi surface in a driven optical lattice : 
Consider a degenerate spin-1/2 Fermi gas trapped in a 
static 3D optical lattice. Our proposed experimental 
setup is illustrated in Fig. 1(a). An additional ID mov¬ 
ing lattice along the x direction is applied to couple the s- 
and p x -orbital bands of the static lattice. The moving lat¬ 
tice is generated by two counter-propagating lasers with a 
frequency difference of to that matches the s-p x band gap, 
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FIG. 1. (Color online) (a) An illustration of the experimental 
proposal: a moving lattice (orange arrows) induces effective 
two-photon Raman couplings between s- and p^-bands of a 
3D static optical lattice (only show 2D here, gray arrows), (b) 
Plot of different coupling strengths fi, a and /3 as a function 
of the static lattice depth Vo- The moving lattice depth V ./ = 
0.8 Er with the recoil energy Er = h 2 /2ma 2 . (c) The single¬ 
particle band structure for Vo = 3.0Er. h = 5.5t 3 and p = 
3 1 3 with t s = 0.111-Er. The colors represent hybrid orbital 
compositions for each momentum state (red for s- and blue 
for p^-orbital states). Dashed lines: bare s- and shifted p x - 
orbital bands without coupling. 


resulting in a two-photon Raman coupling between these 
two bands. All lasers are far-detuned to avoid heating 
from spontaneous emission. The overall time-dependent 
lattice potential can be written as 

V{r,t) = V 0 cos 2 (k L g)+V x cos 2 (k L x+^ , (1) 

r)—x,y,z 

where Vo and V x are the static and moving lattice depths, 
Jcl = 7r /a with the lattice constant a. 

We consider a large static lattice Vo, but a weak mov¬ 
ing lattice V x (i.e., V x <C Vo); therefore only on-site and 
nearest-neighbor tunnelings need be considered and the 
total wave function |\I/} can be expanded in terms of 
the static lattice Wannier functions |\H) = ]TV c S j|s J ) + 
c Pxi\Px)i where j is the site index in the x direction. |s J ) 
and \jP x ) are the s-band and p^-band Wannier functions 
at the j-th lattice site, c S j and c Px j are their annihilation 
operators, respectively. Along the other two directions, 
p-band is not coupled and only s-band is considered and 
their related indices are neglected here for simplicity. 

Under the Wannier basis, we can derive the single¬ 
particle tight-binding Hamiltonian, where the time de¬ 
pendence in the coupling between different orbits could 
be further eliminated using the rotating wave approx¬ 
imation [32], similar to the well-known two level Rabi 
oscillation. The difference from the Rabi oscillation is 
that the two levels here (s and p x bands) have differ¬ 
ent band dispersions. Physically, there are three types 


of possible couplings between s and p x bands, as illus¬ 
trated in Fig. 1(a), with the coupling strengths given 
by U = ■^(s i |sin(2fc i a;)|pj,), a = ^-{s i \cos(2k L x)\p'i+ 1 ), 

and /? = ¥f(s l \sin(2kLx)\p l x ^ 1 ). The first term U de¬ 
notes the coupling of two orbital states at the same site, 
while the last two terms a and /3 are the couplings be¬ 
tween nearest-neighbor sites [33]. The values of f2, a, 
and (3 calculated from the Wannier functions are plotted 
in Fig. 1(b) (see also Fig. SI [32]). /3 is usually small and 
not important for the physics discussed here. 

The resulting time-independent single-particle Hamil¬ 
tonian in the momentum space can be written as 


H 0 (k) 


e s (k) + h n(fc x ) \ 

U(k x ) e p {k) — h J 


( 2 ) 


under the basis (c s (fc), c Px (fe)) T , where n(fc x ) = 

— asin(fc x a) + /3cos(fc x a), e s (fe) = —2t s {cos(k x a) + 
cos(fcj,a) + cos(/c 2 a)] — p and e p (k) = 2t p cos(k x a) — 
2t s [cos{k y a) + cos(fc 2 a)] — p. t s and t p are the nearest- 
neighbor tunneling amplitudes for atoms in the s- and 
p^-orbital states, respectively. 2 h is the energy differ¬ 
ence between to and the band gap A s . p is the chemi¬ 
cal potential. Note that asm(k x a)a x corresponds to the 
band-pseudospin-momentum coupling. 

In the absence of a, Ho(—k) = Ho(k), revealing that 
the single-particle Hamiltonian is symmetric under inver¬ 
sion transformation. This inversion symmetry is broken 
when a and U coexist. A typical single-particle band 
structure, which is asymmetric along fc^-axis, is shown 
in Fig. 1(c). Here we just show the Fermi surface in the 
k y ,z = 0 plane. The Fermi surface is still symmetric along 
k y and k z directions. The orbital and hyperfine-spin de¬ 
grees of freedom of the atoms are independent; there¬ 
fore the coupling between different orbital states does 
not break the spin degeneracy and the hybrid bands are 
spin balanced at any k point. 

Pairing Hamiltonian: Consider a spin-1/2 Fermi gas 
with equal spin populations loaded on such an asym¬ 
metric orbital band. The dominant on-site atom- 
atom interaction between opposite spins can be made 
attractive via Feshbach resonance, similar as regular 
two component Fermi gases [34]. As a good ap¬ 
proximation, the on-site atom-atom interaction can 
take the same form as the time-independent static 
system [32]. In the momentum space, the inter- 
and intra-band interaction term can be written as 
Him = where 

fci + &2 = fc .3 + &4 due to the momentum conservation 
for the two-body scattering process, p and v denote the 
orbital states of two spins. g pv = g J dx\w p (x)\ 2 \w u (x)\ 2 
is the interaction coefficient for two atoms in two orbital 
states (labeled by p and v), and g is the two-body inter¬ 
action strength in free space. To compare the strengths of 
the interactions between two orbital states, we approx¬ 
imate the lattice potential at each site by a harmonic 
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trap, which is a good approximation when the static lat¬ 
tice is not very weak. The relative ratio of the interaction 
strength is found to be g ss : g sp : g ps : g pp = 1 : 0.5 : 0.5 : 
0.75 [32]. Hereafter, we denote g ss = U. 

Under the mean-field approximation, we can rewrite 
the interaction term with the effective pairing between 
atoms. Because the inversion symmetry is broken for 
the single-particle Hamiltonian, the system may favor 
Cooper pairing with a finite center-of-mass momen¬ 
tum between two fermions of opposite spins. For sim¬ 
plicity, the chemical potential is chosen appropriately 
where there is only one simple Fermi surface (see Fig. 
1(c)), therefore we could consider a plane-wave FF-type 
inter- and intra-band pairing A /J1/ (x) = A pu e iC *' x , sim¬ 
ilar as that in spin-orbit coupled system [15-21]. Here 
A pv = g t iu(ofij.(Q/2 + k)c± v (Q/2 — k)) denotes the am¬ 
plitude of the s-wave order parameter between two or¬ 
bital states /i and v, and the FF vector Q = (Q,0,0) 
is the Cooper pairing momentum which is along the 
moving lattice direction. Note that the effective pair¬ 
ing on the asymmetric Fermi surface (Fig. lc) could be 
fc-dependent (i.e., with non-s-wave components) due to 
the fc-dependent hybridization coefficients (determined 
by the eigenfunction of the Hamiltonian (2)) of two 
orbital bands for the asymmetric Fermi surface [35]. 
In the basis of spinor (\H(Q/2 + fc), T*(Q/2 — fc)) T 
with T = (cf s ,cj, s ,c-|- Px , ci Px ) T , the Bogliubov-de Gennes 
(BdG) Hamiltonian can be written as 


HBdG(k) = 


H 0 (y + fe )® cr o 


A* 

^4x4 


^4x4 

~H 0 {Q -k)®a 0 


where <7* (i = x, y, z, 0) are the Pauli matrices, 
A.,, -A, 


A4 X 4 = 


x sp 

A ps App 


( 3 ) 

(~i<7y) . (4) 


For each set of system parameters (Vo, V x , U), the cor¬ 
responding parameters H, a, /?, t s , t p in the BdG Hamil¬ 
tonian (3) are calculated from the Wannier functions, 
from which the order parameter amplitude A pl , and the 
FF vector Q are simultaneously obtained by minimiz¬ 
ing the thermodynamic potential. When A pv ^ 0 and 
<2^0, the system is in an FF phase. When A pv ^ 0, 
Q = 0, the system is in a BCS phase. Otherwise, the 
system is a normal gas. 

Phase diagrams: In Fig. 2 we plot the intra-orbital 
order parameter A ss and the Cooper pairing momentum 
<2 with respect to the static and moving lattice depths 
Vo, V' x . A sp and A pp are much smaller than A ss [32], 
which is ascribed to the initial dominant populations of 
the s-orbital band at the position of the chemical po¬ 
tential. The s- and p^-orbital band tunneling and cou¬ 
pling parameters (t s , t p , H, a) depend on the static lat¬ 
tice depth Vo implicitly; therefore A ss does not change 
monotonically. However, Q and a depend on V' x linearly, 



FIG. 2. (Color online) Phase diagrams of FF superfluids. The 
color describes the amplitude of (a) the order parameter A as 
and (b) the FF vector Q. Other parameters are U = 6.0U, 
g = 10.0t s , h = 8.0 1 3 . 


which directly determine the single-particle band struc¬ 
ture; therefore <2 increases with increasing V' x . Because 
the coupling between s- and p^-orbital states does not de¬ 
pend on spins (the internal states) of atoms, H, a and /3 
modify the energy dispersion in the same way for the two 
spins, leading to the spin degenerate asymmetric Fermi 
surface as shown in Fig. 1(c). Such spin-balanced asym¬ 
metric Fermi surface has little effect on suppressing the 
order parameter, in contrast to the strong suppression of 
the finite momentum pairing order induced by an exter¬ 
nal Zeeman field. Therefore, A ss is large and does not 
change much in the whole parameter region. <2 is pro¬ 
portional to both and a as shown in Fig. 2(b). When 
V x = 0, all the coupling coefficients vanish and the band 
inversion symmetry is preserved; thus the superfluid be¬ 
comes a conventional BCS state. 

When the on-site interaction U is tuned by changing 
the s-wave scattering length through Feshbach resonance, 
the system undergoes a BCS-BEC crossover. BCS-BEC 
crossover physics of Fermi gases has been widely studied 
in free space and in lattices [36, 37]. Here, we present 
the phase diagram in the U-V x plane in Fig. 3. From 
Fig. 3(a), we see A ss is mainly determined by the U and 
changes only slightly with the increase of moving lattice 
depth V' x . FF states with large Q exist in a large param- 



FIG. 3. (Color online) Phase diagrams in the BCS-BEC 
crossover. The color describes the amplitude of (a) the or¬ 
der parameter A ss and (b) the FF vector Q. Vo = 4 Er with 
t s = 0.0855-Er. Other parameters are g = 10.0f s , h = 8.0 t s - 
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eter region which dominates when the moving depth is 
large. In the weak and medium interaction regimes, the 
order parameters are small and the ground state is mainly 
governed by the single-particle Hamiltonian. Therefore Q 
could be significant because of its sensitivity to the single¬ 
particle dispersion. In the very strong interaction regime, 
the fermions form tightly bound molecules and the influ¬ 
ence of the asymmetric energy dispersion on Cooper pairs 
is negligible. Therefore, Q gradually decreases for a large 
U and the ground state eventually becomes a BCS state. 

Stability of FF superfluids: The stability of the FF 
superfluid may be characterized by the thermodynamic 
potential difference Ep f — Ebcs between the FF ground 
state and the possible BCS excited state (by enforc¬ 
ing Q = 0), which is shown in Fig. 4(a). The larger 
|Uff — -EbcsIj the FF state is more stable. When Vf. = 0, 
the inversion symmetry is preserved and the FF super¬ 
fluid becomes the BCS state, therefore Epp = Ebcs- 
With the increasing Vf, Epp — Ebcs becomes negative, 
indicating that the asymmetric energy dispersion favors 
FF superfluids. In Fig. 4(b), we plot the thermodynamic 
potential E in the A SS -Q plane for the premium values 
of A sp and A pp that minimize the total energy, which 
shows that the FF state is indeed the global minimum of 
the thermodynamic potential. 

We also calculate the energy difference of the FF states 
and BCS states in a 3D lattice when only an external 
Zeeman field is applied to break the spin degeneracy. We 
find that the energy difference in our spin-balanced model 
(~ 3 x 10~ 2 f s ) is much larger than that in the spin- 
imbalanced system (~ 10 -3 f s ) with the same interaction 
strength. Moreover, the FF states only exist in a very 
narrow Zeeman held parameter region (~ 10~ 2 f s ) in the 
spin-imbalanced schemes and thus it is hard to find their 
signature experimentally. In contrast, the FF superfluids 
in our spin-balanced system induced by the driven optical 
lattices exist in almost the whole parameter region. 

Q and a are proportional to the moving lattice depth 
Vf, and can be tuned in a wide parameter range to 
achieve an extremely asymmetric energy dispersion. In 
Fig. 4(a), we see Epp — Ebcs decreases sharply when Vf. 
is large. Therefore with a larger Vf, \Eff — Ebcs \ may be 
much larger than that shown in Fig. 4(a), which is gener¬ 
ally impossible in the spin-imbalanced Fermi gases. This 
advantage, together with the large parameter region for 
FF states, make our proposed spin-balanced Fermi gas 
experimentally more feasible for observing FF superflu¬ 
ids than the spin-imbalanced systems. 

Experimental observation: The proposed FF super¬ 
fluids can be realized with different types of fermionic 
atoms, such as 40 K and 6 Li. In the following, we illus¬ 
trate the experimental setup and observation using 40 K. 
The ultracold 40 K gas with a spin-balanced mixture of 
internal states |F,mp) = |9/2, — 9/2) and |9/2, —7/2) 
[52] is trapped in a 3D static optical lattice created by 
counter-propagating far-detuned lasers with wavelength 



FIG. 4. (Color online) (a) The thermodynamic potential 
difference of the FF states and the possible BCS state as 
a function of V/ for Vo = 4 Er (blue solid line, for which 
t 3 = 0.0855 Er) and Vo = 5 Er (orange dashed line, for which 
t 3 = 0.065875 r). (b) The contour plot of the thermodynamic 
potential E in the A SS -Q plane for Vo = 4.0J Er, V/ = 2.0 Er. 
The cross symbol corresponds to the self-consistent solution 
A ss = 3.14fs, Q = 0.377fcz,. Other parameters are U = 8.0fs, 
p = lO.Ofs, h = 8.0t s . 


A = 1064 nrn that defines the wavevector Icl = 2tt/\ 
and the recoil energy Er = H 2 k\/2m = 2 ttH x 4.5 kHz. 
The lowest two orbital bands, s- and p^-orbital, have a 
gap A g ps 2.6 Er when the static lattice depth is tuned 
as Vo = 3.0Er. The ID moving lattice, created by an¬ 
other two counter-propagating lasers with a slight fre¬ 
quency difference of ui ~ A g /h, can be tuned to have a 
lattice depth of Vf = 0.1 ~ 0.8 Er. With these parame¬ 
ters, the resulting coupling strengths have a range of Q 
= 0.24 ~ 1.88f s and a = 0.19 - 1.56f s (t s = 0.111 Er). 
The maximum value of FF momentum Q could be as 
large as 0.3 /cl and the corresponding order parameter 
A ~ t s . For 6 Li atoms, the recoil energy becomes 
Er = 2ttH x 30.0 kHz for the same 1064 nm lasers and 
the relation of all other parameters with Er are the same 
as 40 K. Signatures of FF superfluids can be captured by 
the atom shot noise [38], or the sound speed measure¬ 
ment [39-41], 

Discussion: We note that recently similar moving lat¬ 
tices have been used to couple s-orbital states of neigh¬ 
boring sites of a tilted optical lattice for the simulation of 
the magnetic Hofstadter model with a BEC [42, 43]. In a 
more general context, similar time periodic modulation of 
the lattice to couple different Bloch bands ( e.g ., shaken 
lattices), known as “Floquet engineering” [44-50], has 
been investigated extensively in experiments, leading to 
the observation of various phenomena such as the param¬ 
agnetic to ferromagnetic transition (bosons) [51] and the 
simulation of topological Haldane model (fermions) [52] , 
where the atomic spin states are irrelevant. However, 
the effects of s-wave interaction between two spins of the 
Fermi gas has not been well explored and our proposed 
FF superfluids showcase the rich quantum phases that 
may be generated by the s-wave two-body interactions in 
such Floquet systems. Furthermore, the proposed band- 
pseudospin-momentum coupling in optical lattices may 
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open a new avenue for exploring exotic spin-orbit cou¬ 
pling physics. 
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Note added: Recently, the modification of the band 
structure by the moving lattice has been observed exper¬ 
imentally with a BEC [53]. 
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